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Abstract The path integral of a gauge theory is studied in Coulomb-like gauges. 
The Christ-Lee terms of operator ordering are reproduced within the path integration 
framework. In the presence of fermions, a new operator term, in addition to that of 
Christ-Lee, is discovered. Such kind of terms is found to be instrumental in restoring 
the invariance of the effective Lagrangian under a field dependent gauge transforma- 
tion, which underlies the BRST symmetry. A unitary regularization scheme which 
maintains manifest BRST symmetry and is free from energy divergences is proposed 
for a nonabclian gauge field. 



I. Introduction 



The quantization of a gauge field in the continuum requires gauge fix- 
ing. The standard gauge for perturbative calculations is the covariant gauge, 
in which the theory can be regularized and renormalized systematically. On 
the other hand, noncovariant gauges, though computationally more involved, 
possess a number of advantages of their own. Take Coulomb gauge as an 
example, all degrees of freedom there are physical and unitarity is manifest. 
The Coulomb propagator alone gives rise to the renormalization of the cou- 
pling constant [1] [2]. In addition, the explicit solubility of the Gauss law 
constraint in Coulomb gauge might make it easier to construct variational 
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wave functional to explore the nonperturbative physics of a nonabelian gauge 
field. 

The difficulties with noncovariant gauges of a nonabelian gauge theory 
include 1) the complication of the curvilinear coordinates; 2) the lack of 
manifest BRST and Lorentz invariances; 3) the additional divergence, the 
energy divergence, because of the instantaneous nature of the bare Coulomb 
propagator and the ghost propagator; 4) the lack of a unitary regularization 
scheme which can be applied to all orders and therefore the uncertainty of 
the renormalizability. 

The operator ordering associated with the curvilinear coordinates has 
been investigated by Christ and Lee [3]. They derived the path integral in 
a noncovariant gauge from the Weyl ordering of the corresponding Hamilto- 
nian and found that the effective Lagrangian contains two nonlocal terms, 
referred to as Christ-Lee terms, in addition to the classical Lagrangian and 
the ghost determinant. These terms start to show up at the two loop order 
and their diagrammatic implications have been discussed to that order [4] [5] . 
Here, I approach the noncovariant gauge strictly within the path integral for- 
mulation, along the line of Refs. 6 and 7. Starting with the discrete time 
path integral in the time axial gauge, the change of the integration variables 
to other gauges and the continuous time limit are examined carefully and 
the Christ-Lee terms are reproduced. In the presence of quark fields, a new 
nonlocal term of the effective Lagrangian involving fermion bilinears, which 
starts to show up at the one loop order is discovered. 

According to Feynman [8] , a path integral of a quantum mechanical sys- 
tem is the e — > limit of a multiple integral over canonical coordinates on 
a one dimensional lattice of time slices separated by e. The exponent of 
the weighting factor on each time slice is equal to ie times the classical La- 
grangian only if the canonical coordinates are cartesian [2] . This is the case 
of the time-axial gauge of a gauge theory. Even there, the velocities in the 
Lagrangian is the mean velocity between the neighboring time slices instead 
of instantaneous ones. When transforming the integration variables to curvi- 
linear ones, or to other gauges, e.g. Coulomb gauge, one has to keep track 
of all the contributions to the limit e — > 0, which introduces the Christ-Lee 
operator ordering terms in addition to the classical Lagrangian in terms of 
the new coordinates and the corresponding jacobian. For the same reason, 
the classical Lagrangian with the mean velocities is not exactly invariant un- 
der a gauge transformation which depends on canonical coordinates and the 
variation contributes to the limit e — > 0. The standard form of the BRST 
identity is only recovered after including the Christ-Lee terms. Alternatively, 
one may retain a nonvanishing e and this lends us to a gauge theory with 
discrete time coordinates, which is manifest BRST invariant. For a field the- 
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ory, this discrete time formulation serves as a unitary regularization scheme, 
which regularizes the energy divergence and the ordinary ultraviolet diver- 
gence at one shot. It also possesses several technical advantages which may 
be helpful in higher orders. 

This paper is organized as follows: I will illustrate the technique of the 
gauge fixing within the path integral formulation and the derivation of the 
BRST identity of the soluble model of Ref. [5] in the next two sections. The 
application to a nonabelian gauge field is discussed in the Sections IV and V. 
There I will also test the discrete time regularization scheme by evaluating 
the one loop correction to the Coulomb propagator. The comparison of my 
regularization scheme with others and some comments on the renormalizabil- 
ity will be discussed in the final section. Except for the fermionic operator 
ordering term, the interplay between the BRST identity and the operator 
ordering terms and the discrete time regularization scheme, all other results 
are not new. But it is instructive to see how the operator ordering terms 
come about without referring to the operator formulations. It is also amaz- 
ing to see how similar in formulation the soluble model and the nonabelian 
gauge field are. 

II. The Path Integral of a Soluble Model 

The soluble model proposed by Friedberg, Lee, Pang and the author [5] 
provides a playground for the investigation of gauge fixing and the BRST 
invariance in some nonstandard gauges within the path integral formulation. 

The Lagrangian of the soluble model is 



L = \[{x + giy) 2 + (y - g&f + (i - £) 2 ] - U(x 2 + y 2 ). (2.1) 



It is invariant under the following gauge transformation 



x 



x 



x cos a — y sin a, 



(2.2) 



y^y' 



x sin a + y cos a, 



(2.3) 



z 



z 



1 

z H — a 
9 



(2.4) 



and 
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(2.5) 
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with a an arbitrary function of time. The Lagrangian (2.1) does not contain 
the time derivative of £ and the corresponding equation of motion reads 

dL 

— = g[y(x + g£y) - x(y - g^x)} - i + f = 0, (2.6) 

which is the analog of the Gauss law of a gauge field. In the following, we 
shall review the canonical quantization in the time-axial gauge, i. e., £ = 0, 
convert it into a path integral and transform carefully the path integral into 
the A-gauge, i. e., z = Xx, an analog of the Coulomb gauge. 
In the time-axial gauge, the Lagrangian (2.1) becomes 

L = -(X 2 + Y 2 + Z 2 ) -U(X 2 + Y 2 ). (2.7) 
The canonical momenta corresponding to X, Y, and Z are 

P x = * = - i ^ ( 2 - 8 ) 



and 



Py = Y= -i± (2.9) 
Pz = Z = -i A (2.10) 



and the Hamiltonian operator reads 

H = ^(P x +Py + p z) + U(X 2 + Y 2 ). (2.11) 

The physical states in the Hilbert space are subject to the Gauss law con- 
straint, i.e. 

[P z + g(XP Y -YP x )}\ >=0, (2.12) 

as follows from (2.6) and the operator Pz + g(XP y — YPx) commutes with 
H. In terms of polar coordinates, X = pcos$ and Y = psin$, the wave 
function of a physical state takes the form 

<X,Y,Z\ >=V(p,$-gZ). (2.13) 

For a harmonic oscillator potential, U = \u) 2 {X 2 + Y 2 ), the energy spectrum 
is given by 

E = u(n+ + n- + 1) + -g 2 (n + - n_) 2 (2.14) 

with n+, n- non-negative integers, and the corresponding eigenfunction can 
be expressed in terms of Laguerre polynomials. 



Following Feynman, the transition matrix element < X,Y,Z\e iHt \ > 
can be cast into a path integral 

< X, Y, Z\e- lHt \ >= lim (— ) 1 I J] dX n dY n dZ n 

n=0 



xe-El'o^W <X ,Y ,Z \ >, (2.15) 
where e = t/N and 

L(n) = \{Xl + Yl + Zl) - U{X 2 n + y n 2 ) (2.16) 

with X n = (X n+ i — X n )/e etc.. In the rest part of the paper, the limit sign 
and the normalization factors like (2ine)~ : 2 N will not be displayed explicitly. 

As was pointed out in [3] and [5], the path integral (2.15) picks up the 
velocity 

(X n ,Y n , Z n ) as large as o[e~^j. (2.17) 

In other words, the contribution to the path integral comes from paths which 
can be more zigzag than classical ones. This has to be taken into account in 
variable transformations. The magnitudes of X n , Y n and Z n , on the other 
hand, remains of the order one with a well defined initial wave function 
< Xo,Yq, Zq\ >. To transform the path integral (2.15) to A-gauge, i.e., 
z = \x, we insert the identity 

f N-l 

1 = const. / Y[ d9 n J n S(z n - Xx n ), (2.18) 

J . n 



with 



and 



we have 



71=0 



Xn = X n cos 9 n Y n sin $ n , (2.19) 

y n = X n sin n + Y n cos 9 n , (2.20) 

z n = Z n + -9 n (2.21) 
9 

Jn = ~ + Ay n , (2.22) 
9 



f N-l 

< X,Y,Z\e~ iHt \ >= const. / ] [ dX n dY n dZ n dd n J n 5{z n - Xx n )x 

J n=0 
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xe *E^(») < X ,Y ,Z \ >, (2.23) 
Introducing back £ n via 

in = -k = 9n+1 - 9n (2.24) 
9 9e 

and changing the integration variables from X n , Y n , Z n and 9 n to x n , y n , z n 
and £ n , we obtain that 



< X, Y, Z\e lHt \ >= const. \\_ dx n dy n dz n d£ n 5(z n — Xx n )x 

J n=0 

xe^EnJ 1 '^ < x ,yo,z \ >, (2.25) 



where 



L'(n) = L(n) - -\nj n (2.26) 

with L(n) the same Lagrangian (2.16). Written in terms of the new variables, 
L(n) becomes 

L (n) = ^ (rn+ie"^^ 2 - f n )(e^ 2 r n+ i - r n ) + \{z n - in? ~ U(r n r n ), 

(2.27) 

where we have grouped x n and y n into a 2 x 1 matrix 

rn =(;») (2.28) 

and o"2 is the second Pauli matrix. 

For finite e, (2.25) with (2.26)-(2.27) define a one-dimensional lattice 
gauge field and the limit e — > corresponds to its continuum limit. As 
e -> 0, it follows from (2.17), (2.19)-(2.21) and (2.24) that 

(x n ,y n ,zn,£ n ) = 0(e-v). (2.29) 

Several terms beyond the naive continuum limit have to be kept when ex- 
panding the exponential e te9 ^ n(72 in (2.27) according to £ n [3]. In addition, 
the commutativity between time derivative and path integral contractions 
requires the Lagrangian to be written in terms of x n , y n , z n , £ n , x n , y n and 
z n with x n = (x n + x n+ i)/2 etc.. In another word, we write 

L (n) = L (x n , y n , z n , £ n , x n , y n , z n ) = L (x n , y n , z n , £ n , x n , y n , z n ) + 5L (n) 

(2.30) 



with 

SL (??.) = L (inj Vni Zfii Cni %ni Unj Zn) L (x n) Vm %ni £n? %m Vni z ri)- (2-31) 

According to the estimate (2.29), the contribution from the potential energy 
U to the difference (2.31) vanishes in the limit e — > 0. But that from the 
kinetic energy and from the jacobian (2.22) do not. With this precaution, 
we rewrite L'(n) as 

L'(n) = Leff.(n) + ^(ln J n+1 - \nj n ), (2.32) 

where 

L eS ,(n) = L(n) --lnj n + -\ln J n - hln J n+1 + In J n )] (2.33) 

with 

J n = - + Xy n . (2.34) 
9 

The path integral (2.25) becomes then 

< X, Y, Z\e lHt \ >= const. J N 2 / ] [ dx n dy n dz n d£ n 6(z n - Xx n )x 

J n 



xe^^^-W^ < x , yo^ol >, (2.35) 

Now it is the time to take the limit e — > 0. The small e expansion of L e fj.(n) 
reads 

L eff .(n) = L cL (n)-hnJ n + AL(n) + 0(e5), (2.36) 

where 

^cl.W = ^[(^n + ^n^n) 2 + (yn-^nyn) 2 + (in-^n) 2 ]-^(^n + ^) ( 2 -37) 

is the classical Lagrangian but with mean velocities and 

AL(n) = --e 2 g 2 il{r n -igi n f n a2) [r n + -g<J2£,nf n ) + g e ( 1 + A ^ )2 " ( 2 - 38 ) 

The first term of AL(n) comes from the kinetic energy, the second term from 
the jacobian. It follows from (2.29) that AL(n) = 0(1) and the terms not 
displayed all vanish as e — > 0. 



It remains to convert AL(n) into an equivalent potential (to eliminate 
the explicit e dependence). The recipe was given by Gervais and Jevicki [6], 
which we shall outline here. We assume that the integrations on the time 
slices t = 0, e, 2e, (n — l)e have been carried out and we are left with 

_ 1 f N-l 

<X,Y,Z\e lHt \ >= const. / JJ dx m dy rn dz m d^ m 

m=n+l 

xS(z m - \ Xm )e ie ^^ LcsXm) 

x J dx n dy n dz n dCnJnS(z n -Xx n )e ieLcl i - n) [l+ieAL(n)] < x n , y n , z n \e~ ineH | >, 

(2.39) 

where the corresponding transition matrix element, < x n ,y n , z n \e~ meH \ > is 
a smooth function of x n , y n and z n . The structure of AL(n) is 

AL(n) = J2Cl(x n ,yn)Pi(n)e- (2.40) 
l 

with Pi(n) a product of x n , y n and £ n and n\ the number of factors. Changing 
the integration variables from x n , y n and £ n to x n , y n and £ n while replacing 
(x n , y n ) by (x n+1 - ex n , y n +\ - ey n ) and (x n , y n ) by (x n+ i - ex n /2, y n+1 - 
ey n /2), we have, upon a Taylor expansion in terms of ex n and ey n , that 

(2.39) = const. 2 / jj dx m dy m dz n d£ m 5(z m - Ax m )e* e ^™=«+i cffAmJ 

dx n dy n dz n( % n jh(zn - Xx n )e ieL ^[l - ieV(x n ,y n ) + 0{£)\ 

x <x n ,y n ,z n \e~ ineH \>, (2.41) 

where 

V(x„, y n ) = - < AL(n) > G auss= - ^Cife ?7n) < Pi(n) > Gauss . (2.42) 

and the Gauss average < ... >Gauss is defined to be 

, . f dx n dy n dC n e ieLcl {n) F(n) 

< F( "> *"-= j^^w (2 ' 43) 

while regarding x n and y n constants. Such a procedure is valid even if a 
linear term of x n , y n and £ n with coefficients of the order one is added to 



/ 



L c \ (n), as will be the case when external sources are introduced to generate 
various Green's functions. Introduce a 3 x 1 matrix 




n 

Vn 9 x nin 

in Xx n 



(2.44) 



we have 



1 




1 + ^9Vn 

It follows from (2.43) that 

^ CinCjm ^ Gauss 













lew 



ft- ■ 



(2.45) 



(2.46) 



Working out the Wick contractions in (2.42) according to (2.45) and (2.46), 
we end up with 



V(x,y) 



^(2 + 3A 2 ) + AVy AYV(1 + A 2 ) 



+ 



(2.47) 



8(1 + Xgy) 3 8(1 + Xgy)^ ' 

which agrees with the result obtained via Weyl ordering. The effective La- 
grangian L e g.(n) in (2.35) is then replaced by that of Christ-Lee type, i.e. 



£{n) = L cl (n) - - In J n - V(n). 
e 



(2.48) 



Before closing this section, I would like to remark that the subtleties of 
the path integral depends strongly on the way in which the gauge condition 
is introduced. Consider a general linear gauge fixing with the insertion (2.18) 
replaced by 

N-l 



1 = const 



/n 

J n=0 



d8 n Je- ie ^ {zn ~ Xxn -^ n)2 



(2.49) 



with a and k > gauge parameters like A. This is the discrete version of the 
gauge fixing used in [9] and the A-gauge, (2.18), corresponds to a = and 
k = 0. The analysis in Appendix A gives rise to the estimates in the table I 
for the typical magnitude of £ in the path integral with different choices of 
the gauge parameters. 

For k 0, the limit e — > is trivial and the same estimates apply to the 
gauge fixing with x n and z n in (2.49) replaced by x n and z n , an analog of 
the covariant gauge in a relativistic field theory. But the limit k — > with a 
continuous time will entail higher degrees of energy divergence for individual 
diagrams. 
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Table I. The typical magnitude of £ 





« = 


« ^ 


a = 




0(1) 


a ^ 


o(6-i) 


0(1) 



III. The BRST Identity of the Soluble Model 



There are two approaches to the BRST identity of the soluble model 
(2.1). One can prove BRST invariance by introducing ghost variables and 
establish the identity with external sources. One may also start with the 
Slavnov- Taylor identity [10] and construct the BRST identity afterwards. It 
turns out the former is more straightforward for the path integral (2.25) with 
the lattice Lagrangian (2.26) and (2.27), while the latter is more convenient 
with the Christ-Lee type of path integral. We shall illustrate both approaches 
in the following. 

III.l Prior to the e-expansion 

Introducing the ghost variables c n and c n and an auxiliary field b n , the 
path integral (2.25) can be cast into 

.N-l 

< X, Y, Z\e~ lHT \ >= const. / Y\ dx n dy n dz n d£, n db n dc n dc n x 



xe i£ C=o l brst(«) < XQiy0iZQ \ > (3.1) 
where 



^brstW = ^(r n +ie-^ 2 -^(e^ 

+b n (z n - \x n ) - c n (l + \gy n )cn- (3.2) 

The integration measure and the Lagrangian (3.2) is invariant under the 
following transformation 

5r n = -igO n a 2 r n , (3.3) 

Sz n = 9 n , (3.4) 
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6£n = On, (3.5) 

Sb n = 0, (3.6) 

5c n = (3.7) 

and 

5c n = s n b n (3.8) 

with 9 n = s n c n and s n a Grassmann number. For n-independent s n = s, an 
operator Q such that 5 = sQ can be extracted. It is straightforward to show 
that Q 2 = and the transformation (3.3)- (3.8) is of BRST type. 

To establish the BRST identity, we introduce the generating functional 
of connected Green's functions, 

N 

e tW(JX,u, v ,v) = lim < | e -tffT| >=const . / TT dx n dy n dz n dt n db n dc n dc n x 

T ^°° J n=0 

x < \x N ,y N ,z N > e -EjW(n)+W(n)] < XOiy0iZQ \ ( 3. 9) 
where L ext .(n) stands for the source term, i.e. 

^ext.(n) = J n r n + ( n Z n + U n i n + fj n C n + C n T] ni (3.10) 

| > denotes the ground state of the system, and the limits e — > and 
iVe = T — > oo are understood for the right hand side. It follows from 
the transformations (3.3)-(3.8) that 

< 5L cxt Xn) > f ,= -igJ n (T2 < c n r n > n +(( n - ii n ) < c n > n + < b n > n i] n = 0. 

(3.11) 

This is the prototype of the BRST identity and can be converted into various 
useful forms. 

III. 2 After the e-expansion 

After integrating out the ghost variables and carrying out the e-expansion, 
the path integral (3.9) becomes 

N 



jW(J,t,u,n,fj) = lim < \ e ~iHT\ >=const _ f TT dx n dy n dz n d£ n db n x 

T ^°° J n=0 

x < \x N ,y N ,z N > e «E„[A»)+W(n)] < Xq ^ Zq \ >? (3 . 12) 
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where 

% — 

C(n) = L cL (n) + b n (z n - \x n ) - -In J n - V(n) (3.13) 

is the Lagrangian of Christ-Lee type with L c \ (n) given by (2.37), V(n) by 
(2.47) and L cxt (n) by (3.10) at i] n = fj n = 0. One may introduce the ghost 
variables for the path integral (3.12), but they will be different from the ones 
in (3.9), since the argument of the jacobian J n has been shifted from y n 
to y n . On the other hand, the BRST identity can be constructed from the 
Slavnov- Taylor identity and we shall adapt this strategy. We consider a field 
dependent gauge transformation 

$r n = -iXnWn, (3.14) 

$Zn = ~Xn, (3.15) 

9 

<5£n = ~Xn, (3.16) 
9 

where 

Xn = y~~T (3-17) 
g + Ay n 

with e n an infinitesimal ordinary number. Keeping in mind that the velocities 
%n, Vn and z n , and the coordinates x n , y n and z n follow strictly the discrete 
time definition and the variable transformation (3. 14)- (3. 16) is nonlinear, 
the variation of the Lagrangian L c \ {n) contributes to the path integral in 
the limit e — > 0. We find that 

SL cL (n) = -e g tnXn(Dr) n r n = — e n , 3.18 

4 4(1 + \gy n y 

where the terms containing e n have been dropped since e n is assumed smooth 
with respect to n i.e. e n = 0(1). Furthermore, the combination dx n dy n dz n d^ n db n J n 
ceases to be invariant like the combination dx n dy n dz n d^ n db n J n . We have, 
instead, 



S ( dx n dy n dz n di n ) = dx n dy n dz n d£ n A r 



with 



il la 



(1 + \gy n ) 2 (1 + Xgy n +i) 

1 e\ 2 g 3 x n y n 1 e\ z g A x n yl \ 
+ 4(1 + A^ n )3 2(l + \gy n )*r n - [6AU) 
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With these additional terms, the identity (3.11) is replaced by 



+ < SL c i(n) - iA n - SV(n) > 



= 0, (3.20) 

where the term < b n > comes from the variation of the gauge fixing term 
of (3.13). Upon utilizing (2.45) and (2.46) for the Gauss average in the last 
section, we obtain that 



J^£ n < 8L c i(n) - iA n >=^2e n < - *A 3 # 4 



| 3 A/(2 + 3A 2 )x n -AV(2 + A 2 )x n y n 1^6^ , ^ 



8 (l + Xgy n ) 4 

.3 



8 



(i + Xgyn) 5 ' ' 2 /vy v ^ /v ; i + A^ n ) 6 

= $>»[< SV(n) > +0(ev)], 



> 



+0(e2) 
(3.21) 



and the Slavnov- Taylor identity follows in the limit e — > 
-igJ n CT2 < 



' r n > +(Cn - Un) < 1 , ^ > + < b n >= 0. (3.22) 



(1 + Asry n ) 



i + ^gy n 



This can also be obtained from (3.11) after integrating over c and c at 1] = 
fj = Q. 

To construct the BRST identity, we introduce the ghost variables by 
rewriting (3.12) as 



N 



e iW(J,C,u, v ,fj) = lim < | e -iffT| >= congt _ / | J dx n dy n dz n d^ n db n dc' n dc' n 



T^oo 



X 



n=0 



X < \%N, UN, Z N 



> g^Elo^'W+^xt.H] < ^ ^ > (3.23) 



where 



C\n) = L c \_ (n) + b n (z n - \x n ) - c' n (l + Xgy n )c' n - V(n) (3.24) 



and 



Lext.(n) = J nfn + CnZ n + «n£n + Vnc' n + c' n T] n . 

13 



(3.25) 



Note that the primed ghosts are different from the original ones. Denoting 
the average with respect to the path integral (3.23) by < ... > v , we have, for 
a function of the integration variables, F, 

< F >"= *T(nd«in) ■ (3 - 26) 

Then it follows that 

< c' n r n > v =< 1 _ r n > rjn, , (3.27) 

< c - ><=< rrk: > % (3 ' 28) 

and 

< b n > v =< b n > . (3.29) 

In the limit e — > 0, the difference of y n in (3.27) from y n may be neglected. 
The Slavnov- Taylor identity (3.22) implies the following BRST identity 

-igJ n CT2 < c n r n > v + (Cn - u n ) < c n > v + < b n > v r\ n = 0, (3.30) 

which is equivalent to (3.11). As will be shown in Appendix B, the invariance 
of the Lagrangian under the field dependent transformation (3.14)-(3.17) is 
related to an symmetry of the corresponding Hamiltonian in the A-gauge after 
factoring out the Gauss law constraint. There we shall present a derivation 
of the Slavnov- Taylor identity (3.22) from canonical formulations. 



IV. The Path Integral of an Nonabelian Gauge Field in Coulomb Gauge 

IV. 1 Quantization in the time axial gauge 

The Lagrangian density of a nonabelian gauge theory is 

L = - J d 3 r[±V^V^ + ^74(7^ + m)1;] , (4.1) 



where 
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V l v= ^_^ +gf lmn v rnyn ^ 



with Vl the gauge potential and f lmn the structure constant of the Lie algebra 
of the gauge group. The fermion field ip carries both color and flavor indices 
and the mass matrix m is diagonal with respect to the color indices. The 
gauge covariant derivative is D u = -^- — igT l V^ with T l the generator of the 
gauge group in the representation to which ip belongs. The normalizations of 

flmn and T l giyen by trT l T l' = ljll' and flmnjl'mn = Q^ll' with q 2 the 

second Casmir of the gauge group. The Lagrangian (4.1) is invariant under 
the following gauge transformation 



V, - VI = uV^ + -u— (4.3) 

and 

qjj — > ijj' = U ip (4.4) 

with Vfj, = V^T 1 and u the transformation matrix in the representation of ip. 

The quantization of the gauge field is specified in the time axial gauge 
where Vq = 0. The Lagrangian (4.1) becomes 

/l 1 
d 3 f[-Vjvj - -B)b) + ztft* - ^^jDj + m)tt] , (4.5) 

and the corresponding Hamiltonian reads 

H = J d 3 f [injnj + \b)b) + *t 74 ( 7jjDj + m )*] , (4 . 6 ) 

where the canonical momentum 

□5(0 = Vj(rl = (4.7) 

and Sj(r) = \eijkVjk{r) is the color magnetic field. The Hamiltonian (4.6) 
commutes with the generator of time-independent gauge transformations, Q l 
with 

gl = l(glm Vj _ gflmnyn^m + ( 4 g) 

A physical state in the Hilbert space is subject to the Gauss law constraint, 
i.e. 

Q l \>=0. (4.9) 
The path integral in the time axial gauge can be readily written down 



< V\e-' tHt \ >-- 



con S t./n[™]^E^<»»<^(0, r 1|>, (4.10) 
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where 



[dVdV(M] n = Y[ dV](n, r)dt>(n, r)d^(n, f), 

r,j,l 

L(n) = J <Pf[±Vj(n)Vj(n) - l -B]{n)B]{n) 
+i*(n)74*(n) - ^(n)(^fjDj(n) + m)*(n) 



(4.11) 



with 



VHn) = ^ 



Vj(n+1)-Vj(n) 



J v > _ 



= 0(e-s), 



(4.12) 



(4.13) 



and the initial wave functional satisfies the Gauss law (4.9). The dependence 
of the field amplitudes on r has been suppressed in (4.12). 

IV. 2 Transformation to Coulomb Gauge 

Inserting the following identity into the path integral (4.10), 
1 = const. J Yl du(n, r)J(n)5(V j A l j (n, r)), 

where 



n,r 



Aj(n, r) = u*(n, r)Vj(n, r)u(n, r) H — u'(n, r)Vju(n, r), 

u(n, r) is a representation matrix of the gauge group and 

J(n) = det(-VjPj(n)) 
with £>j m = 5 lm Vj - gf lmn A n y Introducing 

J(n,r)u(n+l,r) = e ie 9 A oM^ 



u 1 



and 



ip(n, r) = u'(n, r)^>(n, r) 
$(n, f) = ty(n, r)u(n + 1, r), 



(4.14) 

(4.15) 
(4.16) 

(4.17) 
(4.18) 

(4.19) 



and transforming the integration variables from Vj(n,r) and u(n,f) into 



A l An, r) and A l (n, r), we obtain that 



< Vie 



-iHt\ 



= const, f JJ[dArf^ / ] n 5(V j 4(n,f))e ie Sn i ' (n) < Aj(0,r)\ >, 

•* n 

(4.20) 
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where 



L'{n) = L{n) - % -\nJ{n) - V(0) J d 3 rlnh(n,r) (4.21) 

with L(n) given by (4.12) and h(n, r) the Haar measure of the integration of 
the group element e legAa<<n ^ with respect to Ao(n,r) 



h(n, r) = 1 - 6 ^A l (n, r)4(n, r) + 0(e 3 g 3 ), 



(4.22) 



which does not have an analog in the soluble model. In terms of the new 
variables, we have 

L(n)= [ <Mtr[£>)£>) - 



1 _ 

-ip(n) [ip(n + 1) — e~ 




with 






£j(n) = ~ 


e iegA (n 


) A j (n+ i)e- ie ^o(n) _ 


and 







l 

— ( 



B j( n ) = o e jki V fc Aj(n) - ViA k {n) - ig[A k (n), A^n)} 



(4.24) 
(4.25) 



The Lagrangian (4.23) with (4.24) and (4.25) defines a gauge theory in a 
spatial continuum and on a temporal lattice. The action e£ n L(n) coin- 
cides with the naive continuum limit of the spatial links of Wilson's lattice 
action. But it comes naturally from the definition of a path integral and the 
procedure of gauge fixing. It follows from (4.13), (4.15) and (4.17) that [3] 



(A t j (n,r),A l (n,r))=O(e-2) 
The path integral (4.20) can be rewritten as 



(4.26). 



< V\e- iHt \ > 



const. J-*(N) / JJ[eL4e2V#'] 

J n 

*e£„Wn)j4(0) < A(0)\ > 



(4.27) 



with 



L off .(n) =L(n)-- [In J(n)+\n h(n)} + - [ In J(n) -hln J(n+ 1) +ln J(n)) 

(4.28) 
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where 



Jin) = det(-VjPj(n)) (4.29) 

with Vf = S^Vj-gf^A^n) and Afa) = ^[A l j (n + 1) + A l j (n)]. The small 
e expansion of L e ff X n ) reads 



Leff.(n) = Lcl.(n) - - In J(n) + AL(n), 



(4.30) 



where 



L cL (n) = y <fr^tr[5j(n)5j(n) -^(n)Bj(n)] 
+iip(n) [74(^(71) + z£rAo(n)"i/>(n)) - (^jDj(n) + m)ip(n)] 



with 

and 
AL 



Bj{n) = -Aj(n) - VjA (n) - ig[A (n), Aj(n)}, 



(4.31) 
(4.32) 



(„) = / 



d 3 f 



\g\ 2 f lml ' f akl ' £]{n)A%{n)A k Q {n) [i» + ±tif (n)A b Q (n) 



^^ 2 (r,Z|[V^<n)]- 1 ri7(n)V / [VA(n)]-H m, i^\n)V i Hr,0 



+ ^e5 3 (0)C 2 /4(n)4(n) + ^^(nj^fnl^nl^fn)] +0(e^) 

(4.33) 

with t l the generator in the adjoint representation, [t l ) ab = if alb . The first 
term of the integrand of AL(n) comes from the e expansion of the color 
electric field (4.24), the second term from the shift of the Jacobian J(ri), 
the last term of (4.28), the third term comes from the Haar measure h and 
the last term from the e expansion of the fermionic part of (4.23). We may 
notice the close resemblance of the first two terms of (4.33) with (2.38). 

IV. 3 Converting AL into an equivalent potential 

Following the recipe of Section III, the potential energy which is equiva- 
lent to AL(n) in the limit e — > is 

V = - < AL(n) > Gauss 
/ Uri i dAUn, 7)dA l Q (n, r)5{V 3 A l An, r))e le I ' f^iM^iM AL(n) 



(4.34) 
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while regarding Aj(n, r) constant. The Gauss average of a product of Aj and 
Aq can be decomposed by Wick's theorem. We have 



< £f (f)^) > Gauss = -<^(r, a^, b) = l^rt 3 (f " (4-35) 



< Ag(r)4(0 > Gauss = --(r,a\GV z G\r*,b), 



(4.36) 



< A%{f)A){?) > Ga 



(f,a|GV,-|r , 6) + (r, alGV^XMf', 6) , (4.37) 



< A1{r)A b A?) > Gauss = - SijiPPi?- 7) + (f, alV^r', b) 



+(r, o|X> i GVj|f' , 6) + (f, a|P i GV 2 G'P i |r', 6) 



(4.38) 



where we have suppressed the n-dependence and G = (—VjVj) 1 with V 
from here on to the end of the section defined at Aj(n,r). Substituting 
(4.35)-(4.38) into (4.33), we obtain 

V = -^C 2 g 2 5 3 (0) J ' d 3 f(f,l\GV 2 G\f,l) 

+ ^ g 2 f kam f nal J ^ l\ GV ^ k) ft m | GVj |f, n) 

-\g 2 f kam f nbl J d 3 f J rfV(f,/|GV i |^,A;)(f,n|V i G'|r y ,m)(f,6|P J G'V l |r y ,a) 
+^9 2 f kam f nbl J d 3 r J rfV(f,/|GVi|r ,k)(r*,m\GV j \r,n){r',a\V i GV 2 GV j \r, 
+^C 2 g 2 5 3 (0) J d 3 r(r,m\GV 2 G\r,m) 

+ ^ g 2 f lka fm na J ^ k \ GV ^ /) ft m \ G W j\f, Tl) 

-\g 2 f nka f lma J rf 3 r(r, k\GVj\r, l)(r, m\GV 'j\r, n) 



I 2 yZmfc ^anfc 
12 



y d¥ (r, 



m|GVj|r, /)(f, a|X>jGV 2 G|f, n) 



+(f, a\VjGVj\f, l)(r, m|GV 2 G|f, n) + (f, n|GVj|r, /)(f, a|P i GV 2 G|f, m) 
1 



~g 2 I d 3 r(r,l\GV 2 G\r,m)ij(r)<y 4 T l T m i;(r). 



(4.39) 
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The first term is the Wick contraction of the Haar measure term of (4.33), 
the second to the fourth terms are from the jacobian of the gauge fixing, i.e., 
the second term of (4.33), the fifth to the eighth terms are from the color 
electric field energy, i.e., the first term of (4.33) and the last term is from the 
fermion part. This lengthy expression can be simplified with the aid of the 
following two Jacobian identities, i.e. 

f abc J d 3 r[(r,a\Vi\X)(r,b\Y)(r,c\Z) + (r,a\X)(r,b\Vj\Y)(r,c\Z) 

+ (r,a\X)(r,b\Y)(r,c\V J \Z)] = 0, (4.40) 

[3] and 

jdbl jckl _|_ jbcl j-akl _|_ jcal jbkl q ^ ^-|^ 

First of all, the seventh term of (4.39) is already of the form of Christ-Lee's 
Vi. The covariant derivative T>j of the third term may be moved into the 
middle factor of the integrand according to (4.40), and the result will cancel 
with the second and the sixth terms through (4.41). Upon repeat applications 
of (4.40) and (4.41), the first, forth, fifth and eighth terms will combine into 
Christ-Lee's V2. We have finally 

V = Vi + V 2 + V 3 , (4.42) 

where 

Vi = ^g 2 J d 3 r{r, l'\GVj\r, l)(r, m\GV/t l \r, m), (4.43) 

V 2 = -\9 2 J d 3 rj d 3 r'(r',l , \(5 i , i + V i ,GV i )\f,n)(f,l\(8 ii/ + P i GV i /)|rW) 

x(r,n\t l GV 2 Gt l '\r*,n') (4.44) 

and 

V 3 = -^g 2 J d 3 f^(r)-f 4 T l T m ij(r)(rJ\G\/ 2 G\f,m). (4.45) 

The terms Vi and V2 are the Christ-Lee operator ordering terms for a pure 
gauge theory. The term V3 is new and its expansion in g reads 

V 3 = ~\g 2 j d 3 r[5 lm (r\V- 2 \r) + Zg 2 f ll ' k f knm {r\V- 2 A^V iV~ 2 A™V jV^r) 



+0(g 3 A 3 ) iJj(r)-f 4 T l T m ij(r), (4.46) 

where the term linear in A l j vanishes because of the transversality. Operator- 
wise, this term stems from the normal ordering of the four fermion coupling 
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in the color Coulomb potential, which is necessary for the passage from the 
canonical formulation to the path integral. The details will be explained in 
Appendix C. The effective Lagrangian in the path integral (4.27) is replaced 
by the following Lagrangian of Christ-Lee type in the limit e — > 

C(n) = L cL (n) - % - \nj(n) - Vi(n) - V 2 (n) - V 3 (n), (4.47) 

The formulation of this section for the Coulomb gauge can be easily 
generalized to an arbitrary noncovariant gauge introduced in [3] 

J rfV(f, Z|r j |f', l')A l -(r') = 0. (4.48) 

Since e is the only dimensional parameter in the formal manipulation of 
this section, it would be expected that Ao(n,r) = 0(e~ 1 ) on dimensional 
grounds, different from the estimate of £ for the soluble model and the esti- 
mate (4.26). On the other hand, the field theory in D = 4 suffers from the 
ultraviolet divergences which have to be regularized in order for the path in- 
tegral to make sense. The validity of the estimate Aq(ti, r) = 0(e~2) as well 
as the Christ-Lee path integral depends on an implicit assumption that there 
is an fixed ultraviolet length, which makes the summation over all physical 
degrees of freedom finite, in the process e — > 0. If e is identified with the 
ultraviolet length as in the discrete time regularization scheme of the next 
section, the e-expansion can nolonger be truncated. 

V. The BRST Identity and the Discrete Time Regularization 

Neglecting fermion couplings, the Lagrangian of a nonabelian gauge field 
with discrete times reads 

L(n) = J d 3 rtr[£>)£>) - Bj(n)Bj(n)] (5.1) 

with £j and Bj given by (4.24) and (4.25). The corresponding path integral 
is 

< V\e~ iHt \ >= const. / Y[[dAdbdcdc] n e ie ^n L ^T(n) < y ^ 

J n 

where 

[dAdbdcdc] n = ]~J dA^n, f)db l (n, f)dc l (n, r)dc l (n, r)h(n, r) (5.3) 

f,fj,,l 

21 



and 

£brst(» = L(n)+ j ' d 3 rb l (n,r)V j A l j (n,r)- j 'd 3 f&(n,f^\yjVj(n,r)] ll 'c l '(n,r). 

(5.4) 

The Lagrangian (5.4) and the integration measure of (5.3) are invariant under 
the following transformation: 

SA l j(n, f) = s n V l f(n, r)/(n, f), (5.5) 
5e iegA (n,f) = igSnC \ n ^ ^l^gA^v) _ igSn+lC \ n + 1 ? ^ e iegA (n,r) T l ^ ^ 

5c l (n, r) = -± Sn gf lab c a (n, r)c b (n, f), (5.7) 

5c?(n,r) = s n b l (n,r) (5-8) 

and 

Sb l (n,r) = 0, (5.9) 

where s n is a Grassmann number. For a n-independent s n , a nilpotent charge 
operator can be extracted and the transformation (5.5)-(5.9) is therefore of 
BRST type. Introducing the generating functional of the connected Green's 
functions via a source term, i.e. 

e ,¥ (« = lim < \e~ mT \ > 

T^oo 

= Const. / Y[[dAdbdcdc] n < \A(N) > e ^Ej L BRST(n)+W(n)] < ^| > 
J n 

(5.10) 

with 

-^ext.(^) = 2 J d^rtr[J^{n,r)A^{n,f)+fi{n,r^c(n,f) + c(n (5-11) 

The invariance of (5.3) and (5.4) under (5.5)-(5.9) implies the following BRST 
identity [11] 

J <i 3 rtr[J(n, r)- < Vc(n,r) > — < VqJq^t) > +igfj(n,r) < c 2 (n, r) > 

+ < b(n, r) > 77(71, r)} = 0, (5.12) 

where 

Vc(n, f) = Vc(n, r) — ig[A(n, r), c(n, r)] (5.13) 
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and 

V J (n,r) = j (n,r) + ig[A (n, ), J (n, r)]. (5.14) 
The transformation law of Aq(tl, r), deduced from (5.6), 

SA (n, r) = -V e(n, r) + ^g\ 2 [A (n, f), [A (n, r),9(n, r)]] + ... (5.15) 

has been utilized, only the first term of which contributes to the limit e — > 0. 
The identity (5.12) can be cast into various useful forms [11]. 

Similar to the case of the soluble model, the BRST identity can also be 
constructed from the Slavnov-Taylor identity of the Christ-Lee type of path 
integral (4.27) with L e g (n) replaced by C{n) of (4.47) in the limit e — > 0. 

Unlike the soluble model, the field theory case suffers from an ultravio- 
let divergence which needs to be regularized and subtracted. Owing to its 
manifest BRST invariance, the discrete time Lagrangian (4.23) with (4.24) 
and (4.25) serves also as a gauge invariant regularization scheme with e a ul- 
traviolet cutoff. There are several additional technical advantages with this 
regularization. 1) The energy integration with a continuum time is regular- 
ized by the summation over the Bloch momentum on the temporal lattice. 
This is particularly important for resolving the ambiguities associated with 
the energy divergence. 2) With fixed Bloch momenta on the temporal lat- 
tice, the integration over spatial momenta is less divergent. There is only 
a finite number of divergent skeletons and these can be handled by the di- 
mensional regularization; 3) For fixed lattice momenta, the integrand of each 
Feynman diagram is a rational function of the spatial momenta and can be 
simplified with the aid of Feynman parametrization; 4) Manifest unitarity is 
maintained throughout the calculation. In what follows, we shall test this 
regularization by an evaluation of the one loop Coulomb propagator in the 
absence of the quark fields. 

The expansion of the Lagrangian (5.3) according to the power of g 2 reads 

L BRST (n) = L cl (n) + J d 3 rb l (n,r)V jA\(n,r) 
-J d 3 rc l (n,r)[V ' jVj{n,r)f ' c 1 ' '(n,r) + R n , (5.16) 

where 

- ^ 2 e 2 /^'r fc/ '(^(n)^(n)4(n)i|(n) 
+^V j A l Q (n)A^(n)A k (n)V j A a Q (n)) + ^(0)C 2 g 2 A l (n)A l (n)] (5.17) 



Rn = d 



73 
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a b c 

Figure 1 The one loop diagrams of the inverse Coulomb propagator. 

where at the order g 2 , only terms of an even number of Aq factor are kept. 
The first term of (5.17) comes from the expansion of e ie9A ° in the color 
electric field and the second from the Haar measure. Both of them have 
been included in AL(n) of (4.33). For the reason we shall explain later, 
the perturbative expansion ought to be performed in Euclidean space, which 
amounts to the substitutions e — > — ie, Aq — > — iA^ and Aj — > i-g^- The 
dressed Coulomb propagator reads 

4(^o, fe) = 



k 2 + a(ko, k) ' 

where the one loop contribution to a(/co, k) is given by the amputated Feyn- 
man diagrams of Fig. 1. plus the contribution of (5.17), i.e. 

a{kQ,k) = — ^ Fig. la + Fig. lb + Fig. lc j + contribution from R n 

(5.18) 

with the relevant Feynman rules given in Fig. 2. A wavy line stands for a 
transverse gluon propagator and contributes a factor 

= -/^fe-^) (5-19) 
/cq + k l v k l J 

with ko = 2 sin | and 9 G (— 7r, tt) a Bloch momentum; a dashed line stands 
for a bare Coulomb propagator and contributes a factor 

, xll' 
5 ll d (k) = —. (5.20) 
k l 

A three point vertex of one Coulomb line and two transverse gluons with 
incoming momenta (6*i, k\), (6*2, fe) and (6*3, k%) is associated with the factor 

-zV^sin^^; (5.21) 
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/, Q l k 1 




n 3 k 3 



e k 



m 2 & 2 



/ Q l k 1 < 



n 9 3 k 3 




a Q l k l b Q 2 k 2 

Figure 2 The relevant ingredients of the diagrams 

a three point vertex of two Coulomb lines and one transverse gluon with 
incoming momenta (6>i , fci), (#2, £2) and (#3, £3) is associated with the factor 



-^ mn (%-^)cos|. 



(5.22) 



A four point vertex of two transverse gluons and two Coulomb lines is asso- 
ciated with the factor 



-g\f la ' a f lh ' h + f^) COS % COS ^Sy . 



2 2 



(5.23) 



With these rules, we have 



fig ,a=I^V|£^ S in^/, (5,4) 
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where 

i = J ■0j i <i ij (e\p)d i Ao'W) 



(2?r) Jo y Jo [(p + kz) 2 +Pqx(1 - z) +PqVZ + k 2 z(l - z)] 4 

(5.25) 

with po = - sin |, p = 7 sm 2" an< ^ 4> = @ — 0, k = p — p the external energy 
and momentum. Similarly, 

Fig. lb =C2S ll 'g 2 f^cos" 9 -II, (5.26) 



where 



11= j ■^ s (p+2k),(p + 2k) j d i] {e\p)d ( 0) 
io io \(v + kz 



and 



(2?r) Jo y [(p + kz) 2 +plx(l- z) + k 2 z(l-z)} 3 

(5.27) 

n7r c$ 2 f d 3 p 



Fig. lc=-C 2 ^/_^cos^/^,( % -) 

= -2Ctf« V f ^ cos 2 ? / (5.28) 
y 7_ 7r 27re 2 7 (27r)3p2 + p2 v ^ 



and 



Contribution of # n = -C 2 <r ^ / — / / 



12j_ 7T 2TveJ (27r) 3 y (p2+p 2 x) 2 

(5.29) 

We shall not expose the details of the evaluation of (5.25)-(5.29), but only 
remark on few key points which lead to the final answer. First of all, the 
p-integrations in (5.27)-(5.29) are all linearly divergent, which upon the re- 
placement 

r d 3 p r d D p 

J w ~ J (6 ' 30) 

give rise to Gamma functions with arguments of the form y + integer, and 
therefore yield finite limits as D — > 3. After the ^-integration, the integrand 
for 9- integration is of the dimension of a momentum. Because of the e of the 
denominator of the Bloch momentum po, the leading divergence as e — > is 
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of the order of e , which reflects the usual quadratic divergence. At k = 0, 
we obtain that 

( 5 - 24 ) = ^2 C ^ ll \ (5-31) 
(5.26) = 0, (5.32) 

( 5 - 28 ) = z^ C 29 2 S U \ (5-33) 
(5.29) = -^C 2 g 2 5 11 '. (5.34) 



and 



37r 2 e 2 

If follows from (5.18) that 

a(k Q ,k) = 0, (5.35) 

which renders the net divergence logarithmic. After some manipulations, we 
obtain that 

-(0.*) = -^/P(lnl-|-|ln2) (5.36) 



and the one loop renormalized Coulomb propagator reads 

Z 



d' (k) = I- (5.37) 



with 

Z = l + -^C 2 g 2 (ln — - — -— In 2V (5.38) 
24tt 2 2y V fee 33 22 /' V ; 

the divergent part of which coincides with the charge renormalization [1] [2] . 
We end this section with two technical remarks: 

1). Euclidean time is adapted for the above one loop calculation. This 
turns out to be necessary for the logarithmically divergent diagrams with the 
integration order we followed. Consider a simple integral with a Minkowski 
momentum p = (po,p) 

/d^p % 
(2tt)4 (p2 + m 2 )2 ( 5 - 39 ) 

with p 2 = fp 1 — Pq. If the Wick rotation is performed before the spatial 
integration, the infinite arc, po = Re 1 ^ with R — > 00, < </> < ^ and 
n < 4> < wm n °f contribute. But if the spatial momentum is integrated 
first as we did, the Wick rotation then will pick up a term from the infinite 
arc. As a result, the renormalization constant will be complex unless we 
start with the Euclidean definition of the diagram. 
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2). It may looks puzzling that the very terms of R n which help to cancel 
the quadratic divergence of the one loop diagrams of the Coulomb propagator 
are actually the same terms which contribute to the Christ-Lee anomalous 
vertices which are expected at two loop level. This paradox is tied to the 
identification e with the ultraviolet cutoff. If an independent ultraviolet 
cutoff is introduced for the integration over p and the limit e — > is taken 
before sending the cutoff to infinity, the contribution of R n to the one loop 
Coulomb propagator, (5.28), will vanish as can be seen easily. 

VI. Concluding Remarks 

In this work, we have carefully traced all the subtleties of gauge fixing 
and variable transformation in a path integral of a gauge model, without re- 
sorting to the operator formalism. For a soluble quantum mechanical model 
in A- gauge and for a nonabelian gauge field in Coulomb gauge, the well 
known operator ordering terms are reproduced exactly. In the presence of 
fermionic degrees of freedom, an additional operator ordering term is discov- 
ered. Because of the intrinsic nonlinearity of a BRST transformation, the 
operator ordering terms are found essential in restoring the simple form of 
the identity associated with this transformation. In the field theory case, a 
manifest BRST invariant and unitary regularization scheme is proposed and 
it does give rise to the correct /3-function at one loop order. 

Though this work does not attempt to prove the renormalizability of a 
nonabelian gauge theory in Coulomb gauge, I do not see any problems in 
applying the discrete time regularization scheme to higher orders. The only 
draw back is that the e-expansion of the temporal lattice Lagrangian can no 
longer be truncated since the ultraviolet cutoff is identified with e. 

Alternatively, one may try to renormalize the theory with a Christ-Lee 
type of path integral. Then one has to face the energy divergence and the 
ambiguities associated with it. The coupling with the ultraviolet divergence 
makes it difficult to organize the cancellation in higher orders. Several sce- 
narios have been proposed but none of them [12] goes smoothly beyond two 
loops. On the other hand, the energy divergence is an artifact of the path 
integral, since it is not there with canonical perturbation methods. In prin- 
ciple, one should be able to organize the energy integral before integrating 
spatial momenta and to reproduce the canonical perturbation series. But 
then no advantages of Feynman diagrams have been taken and the path 
integral seems unnecessary. What we need for the renormalization with a 
Christ-Lee type of path integral is an unambiguous scheme which regular- 
ize the spatial loop integral. The only feasible BRST invariant scheme is a 
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spatial lattice. 

At this point, it is instructive to draw some connections of the path 
integral in the continuum with Wilson's lattice formulation [13]. In the 
absence of quarks, the partition function of Wilson's formulation on a four 
dimensional rectangular lattice reads 

Z= f J] dU tJ e^ Sw[U \ (6.1) 

J <ij> 

where Uij is a gauge group matrix on a nearest neighbor link. The simplest 
choice of the action is 

*M = ? £ E *(i " 5 R ^) + SE tl '( 1 - 5 R ^)] • < 6 ' 2 > 

where £/p = UijU jhUuUu for a plaquette P(ijkl) with the subscript s labeling 
the space-like one and t time-like one, a s , at denote the spatial and temporal 
lattice spacings and d the dimension of t/'s. The lattice Coulomb gauge 
condition can be imposed as in Ref. 14. The discrete time regularization 
scheme presented in Section 5 corresponds to the limit at — > after a s — > 
and any regularization corresponding to Christ-Lee path integral follow from 
the limit a s — > after — > 0. 
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Appendix A 



To estimate the typical contributions of x, y, z and £ to the path integral 
in the limit e — > 0, we may neglect the interaction term and consider the path 
integral with the free Lagrangian only, 

L o( n ) = «[*n + Vl + - £n) 2 - (uJ 2 - i0 + )(xl + y£) ~ ~W - ^ n ~ <n) 2 ], 

CL 



(A.l) 
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where the total time interval T = Ne — > oo with N the number of time 
slices between the time interval T, and the infinitesimal imaginary part of 
10 provides a converging factor of the integral. The last term of (A.l) is 
the gauge fixing term (2.49) with the gauge parameter a. Defining the path 
integral average of an arbitrary function of x n , y n , z n and £ n by 

< F >= jY[n dx ndy n dz n d^ n e le ^n L ^F 
JUndxndy n dz n dCne ie ^ L ^ ' 

we obtain the following expressions for various propagators: 

_ 1 r d6 i e -i{n-m)e 
< x n x m > — < V-nUm >— — / 7; Z 9— — ZTTi (AoJ 

e J 7p 2np*p — u z + z0 + 

< ZnZm >= i /* M ^-^-^fv-fl ,-*'-^ (AA) 



e 2tt (1 + np 2 )(l + Kp* 2 )(p*p -u 2 + i0+) 

, t t f dB X 2 p*p + (1 - ap*p) (p*p - IV 2 ) i(n-m)B (* k\ 

^ m e J_„ 2tt (1 + k P 2 )(1 + K P * 2 )(p*p -u 2 + z0+) 1 } 

1 r de X p*e-^ n - m ^ 
_ 1 r de Ae ~i(^-m)e 

< X n z m >- - J ^ ^ {1 + Kp *2 ){p * p _ u 2 + l0+) 



and 



^ > . = _1 r d6 A V - (k P + ap*) (p*p - u 2 ) l{n _ m)e 
n?m " e7_,27r (l + /^ 2 )(l + /^* 2 )(p*p-^ 2 + z0+) 

(A8) 

-i6_-i 

where p = i e . According to the definition of i n and y n , we have 

1 r de ip*pe- i ( n - m ) e 

< X n X m > = < ynVm >= ~ / 7: Z 2 I -n+ ' ( A9 J 

The squares of the typical magnitude of x n , y n , z n , x n , y n , z n and £ n 
inside the path integral in the limit e — > are of the same order as the 
expectation value of their squares, i.e. the propagators (A.3)-(A.5) at n = m. 
It follows from (A.3) and (A.9) that 

< x\ >=< V 2 n>=^ + 0(e 2 ), (A10) 
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but 



< ij >=< >= 1 + | + OCc 2 ) (All) 

for arbitrary A, k and a. On the other hand, the e — > limit of < z 2 >, 
< i 2 > and < £ 2 > are very delicate and we consider the following situations. 
1) None of k or a vanishes, It follows from (A. 4) and (A. 5) that 



2 A Z 



A k 2A k 

K — a- 



2u(l + KUJ 2 ) 2 y/Jzl 1 + KUJ 2 (l + KUJ 2 ) 2 ! 



+ 0(e 2 ), (A12) 



<£n >= 



X 2 u 



X 2 K 2X 2 K 



K — a- 



+0(e 2 ) (A13) 



2(1 + kuj 2 ) 2 4«V«L 1 + /ou 2 (1 + kuj 2 ) 2 } 

and 

< i 2 >= - + finite terms. (A14) 

The small e expansion of the lattice Lagrangian (2.27) is trivial with such 
a gauge fixing. So is the case when x n and z n in the last term of (A.l) are 
replaced by x n and z n . 

2) k = and a — > before e — > 0. This is the A-gauge in the text. It is 
easy to show, using (A. 3), (A. 4) and (A. 5) that 

A 2 

< z\ >= A 2 < xl >= — , (A15) 

i\ 2 

< i 2 >= A 2 < x\ >= h finite terms (A16) 

and 

< £ >= z ( 1 + A ) + finite terms . (A- 17) 

3) « = but a ^ 0. This corresponds to the "smeared A-gauge". We 
obtain from (A. 4) and (A. 5) that 

< zl >= -%- + finite terms, (A18) 
a X 2 

< i 2 >= -2i-^ + i h finite terms (A19) 

and 

a 1 + A 2 

< £ 2 >= — 2i-g + % h finite terms. (A20) 

The equations (A.12)-(A.20) give the announced estimates in Section II. 
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Appendix B 



The Hamiltonian of the soluble model in the A-gauge is given by [5] 



after solving the Gauss law constraint, where 



1 



M^ = J- 2 (y 2 + -,), 



A 



M xy 1 = Mj = J 2 x(--y), 



yx 



9 



Myy 1 = 



and 



Ay + -) 2 + x 2 (A 2 + l) 
9 J 



1 



J = - + Xy. 

9 



(B.2) 

(B.3) 
(BA) 

(B.5) 



It was pointed out that the Hamiltonian (B.l) commutes with the operator 



K = J 1 (xp y - yp x ), 
i.e. [H, K] = [5]. With U = e ieK , we have 



x £ = UxU 1 = x — ej x y 



and 



y £ = UyU 1 = y + ej l x 
for infinitesimal e. Adding the source term 

h(t) = K (t)K + J x (t)x(t) + J y (t)y(t) 



(B.6) 

(B.7) 
(B.8) 

(B.9) 



to the Hamiltonian (B.l), the Schroedinger equation of the state is given by 

d 



i-Q t \t>= h(t)\t>, 



(5.10) 

where the operators follow the time development generated by H, e.g. 

x {t) = e iHt x(0)e- lHt (B.ll) 
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and 

y (t) = e iHt y(0)e- iHt . (B.12) 

The c-number sources K(t), J x (t) and J y (t) are adiabatically switched on in 
the remote past and are switched off in the remote future. It can be shown 
that 

(jvKJ-v^ =J*KJ-v. (5.13) 

with the subscript W standing for the Weyl ordering. The general solution 
of (B.10) reads 

\t>=U(t,t Q )\t Q > (B.13) 

with 

U(t, t ) = T exp ( - i £ dt'h(t')y (B.U) 

Define the generating functional of the connected Green's functions, W(k, J) 
by 

e ,w M =< \U(oo, -oo)| > (5.15) 

with | > the ground state of the Hamiltonian (B.l), the previously defined 
one, W(J, C, u, 77, fj) in (3.12) aXC,=u = rj = fj = {) corresponds to W(0, J). 
Effecting an infinitesimal transformation (B.7) and (B.8) with e(t) — > at 
t — > ±oo, we have 

Q 

i—\t> £ =h £ (t)\t> £ , (5.16) 

where 

\t> £ =U(t)\t> (5.17) 

and 

de 

h £ {t) = (« - — )K + J^M*) + J y {t)y £ {t). (5.18) 

Consequently, 

|t> £ =£/ £ (£,£ )|*0> £ (5.19) 



with 



£/ e (Mo) = Texp(-i^ dt'h £ (t')y (5.20) 



The invariance of the Hamiltonian 5 and its ground state under the trans- 
formation implies that 

< \U £ (oo,-oo) - t/(oo,-oo)| >= 0, (5.21) 
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which, to the linear power of e gives 

fdt{ 

Joe V 



-^<K> t +e(t) [j x (t) < >* - J y^ < 



gx 



1 + Xgy 



>t 



where the canonical average < ... >t is defined as 

< \T[U(oo,-oo)0(t)]\ > 



<0> t - 



< \U(oo, -oo)| > 



= 0, 
(B.22) 



(B.23) 



Converting < \U(oo,—oo)\ > into the path integral and denoting the path 
integral average by < ... > without the subscript t, we have 



< 



< 



gy 



i + ^gy 



>t=< 



>t=< 



gy(t) 



1 + Xgy(t) 

gx(t) 

1 + Xgy(t) 



>, 



> 



(5.24) 
(5.25) 



and 



< K > t U=o = -g < 



x(t)y(t) - x(t)y(t) + Xg[x 2 (t) + y 2 (t)]x(t) 



> U=o- 



1 + g 2 [x 2 {t)+y 2 {t)] 

(B.26) 

The last equality requires some explanation. In the canonical form, we may 
write 

< K > t \ K=0 = J)| K=0 . (5.27) 

On the other hand, the operator K contains the canonical momenta. Per- 
forming a Legendre transformation of the Hamiltonian H + h [15], the term 
of the corresponding Lagrangian which is linear in k reads 



ng- 



xy — xy + Xg(x 2 + y 2 )x 
1 + g 2 (x 2 + y 2 ) 



(5.28) 



and the equality (B.26) follows from (B.27) with the path integral represen- 
tation of W(«, J). Putting back the £ and z, we find 



< K > t U =0 =< £(t) - z(t) > 
and the identity (B.22) becomes 

r*{| <««>-*(«> > +««[■««> < rrliw > ~ m K 



(5.29) 

gxjt) 

1 + Xgy(t) 
(5.30) 



> 



0, 
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For an arbitrary function e(t), the integration sign may be removed after a 
partial integral and the Slavnov- Taylor identity (3.22) at u = ( = emerges. 
The relation 

<&(*)>=^ <£(*)- *(*)>, (5.31) 
which can be checked explicitly, is utilized in the final step. 



Appendix C 



To make the paper self-contained, we shall go through the path integral 
of fermionic degrees of freedom, following the coherent field treatment of the 
Ref. [7]. Consider a pair of fermion annihilation and creation operators, a, 
and a\ with anticommutator 

{a,ot} = l, (C.l) 

The combination a)aa)a is not zero. But if we replace a and a* by a pair of 
Grassmann numbers z and z, the combination zzzz is always zero. Therefore 
there are ordering ambiguities when transforming the canonical formulation 
for fermionic operators to the path integral. The question is which order goes 
through to the path integral simply through the above replacements a — > z 
and a* — > z. We shall discuss the systematics in the following: 

For a system of M fermionic degrees of freedom, represented by the 
annihilation and creation operators aj and aj with 

{a i ,a j } = (C.2) 

and 

{di,at} = 6ij, (C.3) 

we introduce two set of independent Grassmann numbers, z\, Z2,..., zm and 
z\, Z2,--, zm- We also specify that they anticommute with the a's, a^'s 
and commute with the ket or bra of the ground state in the Hilbert space. 
Furthermore, the following integration rule is imposed 



dzj = / dzj = (C.4) 



and 



J Zidzj = J dziZj = 5ij. (C-5) 



35 



Defining a coherent state by 



and its conjugate by 



\z 1 ,z 2 ,...,ZM>=e^ a ^ j \0> (C.6) 



< Zl ,Z2,...,z M \=<0\e^' Zjaj (C.7) 



with |0 > the ground state. It follows that 

0>j\z\, Z2, ZM >= Zj\zi, Z2, zm >, (C-8) 



< Z\, Z2, ZM\a\ =< Zi, Z2, ...Zm\Zj (C-9) 



and 



< z\,z,2, ...,z M \zi,Z2, ...,z M >= e 2 - 1 ^^ 3 . (CIO) 
Furthermore, we have the completeness relation 

/ \z h Z2,...,z M > Y[dzjdzje ZjZj < zi, z 2 , z M \ = 1- (C-H) 
J j 

Let the Hamiltonian of the system be 

H{a\a) = J2 u ij4 a j + \ Yl r "':i'j"i";'"j'".i + - 

ij ii'jj' 



«i,-,'m;3m>-ji 



where the normal ordering with respect to the state |0 > is the crucial point. 
It follows from (C.8) and (C.9) that 

< zi, zm\H\zi, zm >= H(z, z) 

= ^ ^ UJjjZjZj + — ^ ^ fjj' j/jZiZ^ZjiZj + ... 

"*"Xfl Zh—ZiM Z 3M— z 3i (C.13) 

and therefore 



M! 



< zi, ...,^ M |e^|zi, ...,z M >= e^* jZj [l - ieH(z,z) + 0(e 2 )]. (C.14) 
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With the aid of the completeness relation (C.ll), we end up with the follow- 
ing path integral representation of the fermionic system 



. N-l 

< z[, ...,zf M \e- itH \ Zl , ...,z M >= / [dz] N l[ [dzdz] n [dz} e l ^n L ^ (C.15) 



where t = Ne and e — > at fixed t, and we have made the abbreviation 

\dzdz\ n = Y[ dzj(n)dzj(n), (C.16) 
j 

[dz] N = l[dz j (N) (C.17) 

j 

[dz] o = l[d Zj (0) (C.18) 



and 



The Lagrangian L(n) reads 

L{n) = iJ2zj(n)zj(n) - H(z(n), z(n)) (C.19) 

3 

with Zj(n) = \[ Zj {n + 1) - Zj (n)}. 

Like bosonic operators, the ordering ambiguity here is also reflected in 
the difference between the Dyson- Wick contraction and the path integral 
contraction at equal time. Consider a free system whose Hamiltonian is 
given by (C.12) with ujij = u5ij and all u's vanishing. The Dyson-Wick 
contraction gives. 

lim < 0|T(a(t)a t (0))|0 >= 1 (C.20) 

while 

lim < 0\T(a(t)a\0))\0 >= 0. (C.21) 

t->o- 

The path integral, on the other hand, gives rise to an unambiguous result at 
t = since 



_ / [dz] N Yln=i [dzdz] n [dz] QZi(m) zj (m)e te ^« L ^ 



f[dz} N Un^i[dzdz] n [dz] e^ n L ^ 

e J_ n 2n z e - 1 _ w + z o+ 
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To illustrate the caution which is needed in transforming the canoni- 
cal formulation to path integral, we consider a soluble gauge model whose 
Lagrangian is given by 

L = \{z - 2 + i^H^P ~ i9t) ~ mfy (C.23) 
with t/j, ^ fermionic and z, £ bosonic. The gauge transformation reads 

z -> z' = z + -, (C.24) 



rv 

e-e' = e+- (c.25) 

and 

</, -> <// = (C26) 

with a an arbitrary function of time. In the time axial gauge where £ = 0, 
The Hamiltonian corresponding to (C.23) is 

H=-^ + m&9 (C.27) 
and the Gauss law constraint is 

(- < ^-^ t *)l >=0 ( a28 ) 

The constraint can be solved explicitly and the physical spectrum consists 
of two states with ty'ty = 0, 1 and the corresponding eigenvalue of H = 

2 

and m + Though trivial, we still follow the transformation of (C.27) and 
(C.28) to the gauge where z = 0, with the dynamical variables 9 determined 



(C.28) becomes 



by Z + I = and -0 = e****. The Hamiltonian (C.27) and the constraint 



2 d# 2 

and 



2 q2 

^ = -77-579+^^ (C29) 



(-^-VV)|>=0. (C.30) 



Substituting the solution of (C.30) into (C.29), we obtain 



H eS , =m^+ y(^) 2 - (C31) 
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Following the above recipe, we convert (C.28) into a path integral 

< ip\e- itH <*-\ >= I \\dz n d^ n d^ n di) n 5{z n )e ie ^n L ^^ < V> | >, (C.32) 

J n 

where 

L e ftX n ) = 2^ n ~ ^ + i ^ n< ^ n ~ WZn^n) ~ mipnlpn - —^ n ^n (C33) 

with the last term comes from the normal ordering of the four-fermion term 
of (C.31). The integration over £ n in (C.33) will not generate quartic terms 
since the combination ('i/vAn) 2 vanishes. Applying the Feynman rules given 
by the path integral (C.32) at t — > oo, we have verified explicitly that the 
shift of the self-energy because of the interaction vanishes to one loop order, 
in agreement with the result of canonical quantization. 

Finally, we come to the nonabelian gauge field. The four fermion Coulomb 
interaction term of Christ-Lee Hamiltonian in Coulomb gauge reads 

#Coui. = y J d 3 fd 3 f^ij^^T l ij(r)(rJ\G(-V 2 )G\r'j / )ij\r' )T 1 ' ^{r 1 ) 

2 r- 

= 9 — / rfVrfV : iP^r)T l ij(r)(fJ\G(-V 2 )G\^J , )^(f')T l 'ij(r') : 

2 j- 

+y / d 3 r{r,l\G{-V 2 )G\? ' ,l')^{r)T l T l V(r). (CM) 

The last term becomes V3 of (4.45). 

The additional term steming from the normal ordering of fermionic op- 
erators begins to show up at one loop level, unlike its bosonic counterpart. 
In the case of an abelian gauge theory, the term (4.45) corresponds to the 
self Coulomb energy of a fermion and is not observable, but here, for the 
nonabelian case, it carries the coupling to the gluon fields and may not be 
ignored. 
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